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Exact numerical diagonalization is carried out for the Bohr Hamihonian with a /3-soft, axially 
stabihzed potentiaL Wave function and observable properties are found to be dominated by strong /3- 
7 coupling effects. The validity of the approximate separation of variables introduced with the X(5) 
model, extensively applied in recent analyses of axially stabilized transitional nuclei, is examined, 
and the reasons for its breakdown are analyzed. 
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I. INTRODUCTION 

The Bohr Hamiltonian [1] with a /3-soft but 7- 
stabilized potential has served as the basis of recent inves- 
tigations of the collective structure of transitional nuclei 
intermediate between spherical and axially symmetric de- 
formed shape [2]. The solutions obtained so far for the 
X(5) model [2, 3] and its various extensions [4-9] have 
relied upon an approximate separation of variables, in- 
troduced in Ref. [2], since solution of the exact problem 
was not possible. These approximate calculations have 
been extensively compared with experimental data [10- 
20]. 

However, numerical methods recently developed by 
Rowe et al. [21-23] make exact numerical diagonaliza- 
tion of the Bohr Hamiltonian feasible for transitional and 
deformed situations, without recourse to the approxima- 
tions of Ref. [2] . The solution process involves diagonal- 
ization in a basis constructed from products of optimally 



J 



chosen [3 wave functions with five-dimensional spherical 
harmonics. This method yields much more rapid conver- 
gence in the presence of significant (3 deformation than 
is obtained with conventional methods [24, 25]. 

In the present work, an exact numerical solution for 
the X(5) Hamiltonian is obtained (Sec. II). The results 
for wave functions and observables are examined, and the 
properties of /3-soft nuclei are found to be dominated by 
a strong coupling between the /3 and 7 degrees of freedom 
(Sec. HI). Since extensive prior work has been carried out 
using the approximate separation of variables of Ref. [2] , 
this approximation is reviewed and the reasons for its 
breakdown are analyzed (Sec. IV). 



II. HAMILTONIAN AND SOLUTION METHOD 

The Bohr Hamiltonian, in terms of the quadrupole de- 
formation variables /3 and 7 and the Euler angles d, is [1] 
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where L'^ are the intrinsic frame angular momentum com- 
ponents. For the transition between spherical and axially 
symmetric deformed structure, it is appropriate to con- 
sider a potential V{j3, 7) which is soft with respect to 
(3 but which provides confinement about 7 = 0. In the 
X(5) model, a schematic form V{[3,"f) = Vplfi) -\- 1^(7) 
is used, where Vfj is taken to be a square well potential 
[V(3(/3) = for /3 < /3^ and 00 otherwise] and provides 
stabilization around 7 = 6. 

Under the approximate separation of variables of 
Ref. [2], most results are independent of the specific 
choice of V^. Consequently, in Ref. [2], the confining 
potential ^^(7) is simply described as oc 7^ for small 7. 
For solution of the full problem, V^(7) must be defined 



more completely. In general for the Bohr Hamiltonian, 
the potential energy V{(3,^) must be periodic in 7, with 
period 27r/3, and reflection symmetric about 7r/3, to en- 
sure that the potential energy is invariant under relabel- 
ing of the intrinsic axes [25]. The natural choice of such 
potential is Ky(7) oc (1 — cos 37) [Fig. 1(a)], as consid- 
ered in Ref. [21]. However, for consistency with Ref. [2], 
in the present work the 7 potential is chosen to be the 
oscillator potential, given by V^(7) — A-y^ on the interval 
< 7 < 7r/3 and obtained outside this interval from the 
symmetry requirements on 7 [Fig. 1(b)]. 

As usual for eigenproblems involving the Bohr Hamil- 
tonian, the parameter dependence of the solution can be 
simplified by an appropriate choice of dimensionless pa- 
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FIG. 1: Simple forms of the potential V^(7) satisfying the 
coordinate symmetry constraints of Ref. [25]: (a) ^7(7) = 
X(l - cos37) [21] and (b) VCy(7) = (0 < 7 < 7r/3), used 
here for consistency with Ref. [2]. 

rameters. Transformation of the potential as V^(/3, 7) 
c'^V{cP, 7) or multiplication of the Hamiltonian by a con- 
stant factor both leave the solution invariant, to within 
an overall scale factor on the eigenvalues and an overall 
dilation of all wave functions {e.g., Ref. [26]). Conse- 
quently, all energy ratios and transition matrix elements 
obtained from diagonalization of the Hamiltonian (1) de- 
pend only upon the parameter combination 



which measures the "7-stiffness" of the Hamiltonian. In 
the remainder of this article, the notation is simplfied by 

setting H'^/{2B) = 1 and /S^ = 1, so that a = A. (Results 
for any other values of these parameters may then be 
obtained according to simple scaling relations [26].) 

For a = 0, the potential is 7-indcpcndent, and the 
Hamiltonian (1) reduces to the E(5) Hamiltonian [27]. In 
this case, an exact separation of variables occurs [28] , and 
the eigenproblem can be solved analytically [27] . Nonzero 
values of a yield confinement around 7 = 0. The realistic 
range of values for a is discussed in Sec. III. 

Rowe et al. [21-23] have recently proposed a method 
for numerical diagonalization of the Bohr Hamiltonian, 
with respect to an optimized product basis constructed 
from the five-dimensional spherical harmonics [29, 30]. 
Much as a basis for solution of the Schrodinger prob- 
lem in three dimensions can be formed from the prod- 
ucts of a complete set of radial functions fi{r) with the 
three-dimensional spherical harmonics Y^{9,ip), a basis 
for solution of the Bohr problem can be constructed from 
products of radial basis functions fi{f3) with the five- 
dimensional spherical harmonics ^',,aLj\f (7, i?)- 

This method is especially suitable for application to 
transitional and deformed nuclei, since the radial basis 
functions can be chosen to match the particular radial 
potential at hand. It provides vastly more rapid conver- 
gence in the presence of significant (3 deformation than 
is obtained with conventional methods based on diago- 
nalization in an oscillator basis (see Ref. [23]). Also, it 



can be applied to potentials for which the oscillator basis 
methods are simply inapplicable. For the X(5) problem, 
the wave functions must vanish for P> 1. This boundary 
condition cannot be satisfied in a finite basis of oscilla- 
tor eigenfunctions. Instead, in the present work suitable 
basis functions are defined as fi{f3) oc J,y(x^_i/3) for 
/3 < 1 and zero elsewhere, where is the ith zero of 
Ji, . The calculation of Hamiltonian matrix elements with 
respect to this radial basis is described further in the Ap- 
pendix. The choice i/ = 3/2 makes these /»(/?) the exact 
radial wave functions for the seniority zero states in the 
E(5) limit (a = 0) [27]. 

Rowe et al. [22] provide an algorithm for the 
explicit construction of the five-dimensional spheri- 
cal harmonics, as sums of the form 'i'vaLAii"/,'^) = 
EK^oF^^LKil^iiKi^), where 4>MKi^) = ii^L + 

even 

l)/{16nHl+SK))mDi,^{d) + {-)LDi,_^m*. The 
spherical harmonics are seniority eigenstates, labeled by 
the seniority quantum number v, a multiplicity index a, 
and the angular momentum quantum numbers L and M. 
They arc the exact angular wave functions in the 7-soft 
(a = 0) limit of the present problem. In the Hamilto- 
nian (1), the angular kinetic energy operator (the quan- 
tity in parentheses) is simply the negative of the seniority 
operator A [28] . Thus, calculation of its matrix elements 
between spherical harmonics is trivial (A6). The matrix 
elements of an arbitrary spherical tensor function of 7 
and the Euler angles can be obtained through a series 
of straightforward integrations, as detailed further in the 
Appendix. 

High-seniority spherical harmonics are needed for the 
construction of highly 7- localized wave functions. Thus, 
in general, diagonalization for larger 7-stifFncsscs requires 
larger angular bases. For the range of 7 stiffnesses con- 
sidered in the present work (0 < a < 1000), a product 
basis constructed from the first ~ 5 radial functions and 
the 15 to 20 lowest seniority spherical harmonics suffices 
to provide convergence of the calculated observables for 
low-lying states. 



III. RESULTS 

Results obtained from numerical diagonalization of the 

X(5) Hamiltonian, for a = 200 and a = 1000, are shown 
in Fig. 2. Level energies and E2 transition strengths 
for the lowest-lying levels are indicated. More detailed 
tabulations of energy and B{E2) observables, calculated 
for < a < 1000, are provided through the Electronic 
Physics Auxiliary Publication Service (EPAPS) [31]. For 
the transitional rare earth nuclei (with N « 90) the 7 
band head energy is ~ 6 to 8 times the 2^ energy, consis- 
tent with a 7 stiffness a w 200 to 300. This may therefore 
be considered the most "realistic" range of values for the 
7 stiffness parameter, of the greatest phenomenological 
interest. 

First let us note some basic properties of the solu- 
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FIG. 2: Level schemes for a = 200 (top) and a — 1000 (bottom), from the exact numerical solution for the X(5) Hamiltonian. 
Excitation energies of the lowest members of the ground, /3, and 7 bands, normalized to E(2f), are shown at left. Electric 
quadrupole transition strengths, normalized to B{E2;2f — > 0^) = 100, are shown at right. 



tion [Fig. 2]. The levels are organized into clearly- 
defined bands, characterized by strong intraband tran- 
sition strengths. The lowest energy bands have the spin 
contents associated with a rotational ground state band, 
(3 band, and 7 band. Prom a comparison of the a = 200 
and a = 1000 cases shown in Fig. 2, it is apparent that 
such properties as the band head excitation energies, en- 
ergy spacings within the bands, and transition strengths 
are strongly dependent on the 7 stiffness. This is true for 
both the f3 and 7 bands. While it is natural to expect 
the 7 band energy to increase with 7 stiffness, it is seen 
from Fig. 2 that the P band head energy also increases 
with 7 stiffness by a comparable amount. For larger 7 
stiffness [Fig. 2(c)], the energy spacings within the bands 
tend towards rigid rotor [L{L + 1)] spacings. 

Before considering the spectroscopic observables in 
greater detail, let us first examine the underlying fea- 
tures of the wave functions </?(/?, 7, In the approxi- 



mate treatments of the problem [2], the wave functions 
were separable into products of /3, 7, and Euler angle 
functions, which were considered separately. This is no 
longer possible for the full solutions, but we can still 
consider the probability distributions with respect to one 
or more of these coordinates, constructed by integration 
over all remaining coordinates (see Appendix). The in- 
tegration metric appropriate to the Bohr coordinates is 
/3^d/3|sin 3^\d'yd'd [1, 25] . Contour plots of the probability 
distribution with respect to /3 and 7 are shown in Fig. 3, 
for the ground, /?, and 7 band head states. The prob- 
ability distributions with respect to /? or 7 individually, 
i.e., integrated over the other coordinate, are shown for 
these same states in Fig. 4. 

The 7 dependences of the probability distributions of 
the various states [Figs. 3 and 4 (right)] are modulated 
by the |sin37| factor in the integration metric, which 
causes the probability density to vanish at 7 = and 
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FIG. 3: Contour plots of the probability distributions with respect to (3 and 7, from the exact numerical solution for the X(5) 
Hamiltonian, for (top) the ground state (Oj"), (middle) the /3 band head (Oj), and (bottom) the 7 band head (2+). These 

are shown for 7 stiffness parameter values a = Q (left), a = 200 (middle), and a — 1000 (right). Contours of the potential 
[VCy(7) = aj'^] are also shown (dashed curves). For the 7-soft case a = 0, where the bands are not well defined, the 0^ (^ = 1, 
v = 0), O2 (C = 2, « = 0), and 2j (^ = 1, « = 2) states, respectively, are shown. Plots are in the standard polar form, with (3 as 
the radial coordinate and 7 as the angulax coordinate. 



7r/3. In the 7-soft limit, the probability distributions for 
several of the lowest-lying states are simply P(7) oc sin 87 
[Fig. 4(b)] and thus peaked at 7 = 7r/6, or 30°. For the 
7-stabilized cases [Fig. 4(d,f)], the 7 distributions for the 
of and 0^ states are similar to each other, while that for 
the 2+ state is peaked at nearly twice the 7 value. The 
probability distributions arc compressed closer to 7 = 
with increasing Hamiltonian 7-stifFness, but even for the 
stifTest case considered (a = 1000) the value of (7) is ~ 11° 
for the ground state and ~ 16° for the 7 band head. 

Thus, for realistic 7 stiffnesses, the wave functions ex- 
hibit considerable "dynamical" 7 softness. This is per- 
haps contrary to the common conception that nuclei with 
well-defined 7-bands, like those in Fig. 2, are "axially 
symmetric" and have 7 ~ 0. (Recent empirical esti- 
mates of the effective 7 deformation of transitional and 
deformed rare earth nuclei [32] yield comparably large 
dynamical 7 softness.) 



The (3 dependences of the probability distributions of 
the various states [Figs. 3 and 4 (left)] are seen to be 
strongly influenced by f3-j interaction, migrating towards 
larger /3 as the 7 stiffness increases. This evolution can be 
understood in a qualitative fashion in terms of the five- 
dimensional analogue of the centrifugal effect. In the 7- 
soft limit, where the problem is separable, the /3 wave 
function is governed by a radial Schrodingcr equation 
containing a term Vc{l3) = {v+l){v+2)/ [28, 33], analo- 
gous to the centrifugal potential in the three-dimensional 
central force problem, as shown in Fig. 5(a). The 
strength of this term depends upon the seniority quan- 
tum number v, alternatively denoted r in Ref. [27]. The 
1//?^ dependence energetically penalizes small /3 values 
and tends to "push" the wave function towards larger 
,9. But even in the 7-stabilized situation, an analogue 
of the centrifugal effect occurs, and its general features 
may be understood by observing that the angular kinetic 
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FIG. 4: Probability distributions with respect to (5 (left) and 7 (right) for the ground state (solid), the /3 band head (dashed), 
and the 7 band head (dotted), shown for 7 stiffness parameter values (top) a = 0, (middle) a — 200, and (bottom) a — 1000. 
The 7-confining potential Vy and the level eigenvalues are shown in the insets. In panel (b), as in Fig. 3, the 0^ (^ = 1, w = 0), 
Oj (^ = 2, v = 0), and 2j (^ = 1, v = 2) states are shown. These states have identical 7 distributions, ^-"(7) oc sin 87. 
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FIG. 5: (a) Centrifugal potentials K(/3) = {v + l){v + 2)1 f 
for the radial Schrodinger equation in the 7-soft limit, for 
quantum numbers w = and 2. (b) Effective centrifugal po- 
tentials Vc(/3) ~ {{v + l){v + 2))//3^, schematically indicating 
the strength of the centrifugal effect for the 7-stabilized case 
a = 1000. The level energies and centrifugal potentials are 
shown for the same states as in Figs. 3 and 4, namely, the 
ground state (solid), the /? band head (dashed), and the 7 
band head (dotted). 



energy operator [the factor in parentheses in the Hamil- 
tonian (1)] multiplies a /^-dependent factor l//3^. The 
angular kinetic energy contains contributions from both 
the rotational degrees of freedom and the 7 degree of 
freedom. As the 7 stiffness of the potential increases, the 



confinement about 7 = results in a larger 7 contribution 
to kinetic energy, even for the ground state, where this 
is essentially the 7 vibrational zero point kinetic energy. 
Thus, as the 7 stiffness increases, the influence of the 
1 term in the Hamiltonian becomes larger, displacing 
the wave functions towards larger /?. 

An idea of the strength of the centrifugal effect in the 
7-stabilized cases can be obtained by numerically evalu- 
ating the expectation value of the angular kinetic energy 
for any given state. The resulting "effective" centrifugal 
potential V;(/3) = {{v+l){v+2))//3'^ = {k+2)/0^ is shown 
in Fig. 5(b). (This potential is not of calculational value, 
but it is useful in understanding the solution properties.) 
For the square well /3 potential used in the present model, 
the extent of the wave function in f} is limited by the hard 
wall at /3= 1. The centrifugal effect compresses the wave 
function against the wall, so for large 7 stiffnesses the 
wave function becomes localized with respect to /?, just 
within the wall. Thus, even though the potential is flat 
in /3, the /3-7 interaction inherent in the kinetic energy 
operator induces something akin to rigid (3 deformation. 

The 0^ state has a bimodal probability distribution 
with respect to /3 [Fig. 4 (left)]. In the 7-soft limit 
[Fig. 4(a)], an exact zero in the probability distribution 
arises from the node in the radial wave function, but for 
a > an exact zero is not expected. The 0^ state also 
has a smaller mean value of /? than the ground state. 
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FIG. 6: Dependence of energy observables upon the 7 stiffness 

parameter a, for < a < 1000. (a) Excitation energies of 
the (solid), 2 J and 2^ (dashed), and 3]^ (dotted) levels, 
normalized to E{2'^). For large a. those states become the 
lowest spin members of the /3 and 7 bands, (b) Energy ratio 
R^,^ = E{A+)/E{2+). 

The 0^ state is subject to a similar effective centrifu- 
gal potential to that for the ground state [Fig. 5(b)], 
but it has more energy available in the /3 degree of free- 
dom and is thus less strongly confined against the wall 
at /3 = 1. The 2+ state instead has a larger mean (3 
than the ground state. The angular kinetic energy for 
the 7 band is about twice that for the ground state, in 
the limit of harmonic 7 oscillations [2], and hence so is 
the centrifugal effect [Fig. 5(b)]. For a = 200 the mean /? 
values are (/3)o+ w 0.64, k. 0.54, and {I3)^+ k. 0.70, 

while for a = 1000 they are w 0.71, {(3)^+ k. 0.60, 

and (/3)2_+ wO.76. 

The dependences of a few basic observables upon the 
7 stiffness are shown in Figs. 6-8. Naturally, the 7 band 

energy increases with 7 stiffness: the evolution proceeds 
from the E(5) limit (a = 0), through 7 excitation ener- 
gies appropriate to rare earth transitional nuclei (a w 200 
to 300), to very 7-stiff structure [Fig. 6(a)]. But the 
/? excitation energy increases with 7 stiffness as well, at 
about half the rate at which the 7 excitation energy docs. 
An avoided crossing of the 2^ and 2+ states occurs for 
OW450. 

The angular momentum dependence of energies within 
the yrast band varies substantially with 7 stiffness 
[Fig. 8(a)]. For large a, it approaches the L{L + 1) depen- 
dence of the axially symmetric rigid rotor. In particular, 
the energy ratio i?4/2 = -E(4f )/£^(2f ) ranges from 2.20 in 
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FIG. 7: Dependence of B{E2) observables upon the 7 stiffness 
parameter a, for < a < 1000. (a) Branches from the 2j level 
to the yrast 0^ (solid), 2'^ (dashed), and 4'^ (dotted) levels, 
normalized to BE{2f Of)- Immediately below the avoided 
crossing at a w 450 this is the 2+ level, while immediately 
above it is the 2 J level (as indicated), (b) Branches from 
the 2+ level to the yrast 0+ (sohd), 2+ (dashed), and 4+ 
(dotted) levels, normalized to BE{2'^ — > 0^). Immediately 
below the avoided crossing at a w 450 this is the 2^ level, 
while immediately above it is the 2+ level (as indicated). 



the E(5) limit to > 3.2 for large a [Fig. 6(b)]. The yrast 
band B{E2) strengths [Fig. 8(b)] also vary with stiffness, 
for large a likewise approaching rigid rotor values. The 
observables are thus seen to clearly reflect the rigid (3 
deformation induced by the centrifugal effect at large 7 
stiffness. 

Further inspection of Fig. 2 provides a more detailed 
view of the spectroscopic properties obtained. Strong in- 
terband E2 transitions are predicted, some comparable 
in strength to in-band transitions. A radical suppres- 
sion of the spin- descending 13 band to ground state band 
transitions {e.g., 2^ ^ Oj") relative to the spm- ascending 
transitions {e.g., 2'^^ ^ A'l) is one of the most notable 
features. (Only the branching of the 21" state is shown 
in Fig. 2, but the branchings of the higher-spin band 
mcnnbers are provided through the EPAPS [31].) The 
strengths of the interband transitions depend in detail 
upon the 7 stiffness, as shown in Fig. 7. As the 7 stiff- 
ness increases, the [3 to ground band and 7 to ground 
band transitions become weaker overall. They also tend 
closer to the rigid rotor Alaga rule branching ratios [1]. 
Strong /3-7 interband transition strengths occur as well 
[Fig. 2], but these vary greatly with 7 stiffness. 
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FIG. 8: Yrast band (a) energies and (b) B{E2) strengths as 
functions of angular momentum, for a — 0, 200, and 1000 
(dashed). The values obtained under the approximate sepa- 
ration of variables (Sec. IV) are indicated for reference, as are 
the rigid rotor values (solid). 



For all 7 stiffnesses, the energy spacing scale of the 
(3 band is enlarged relative to that of the ground state 
band [Figs. 2 and 6(a)]. (In contrast, the spacing scale of 
levels within the 7 band is similar to that of the ground 
state band.) The enlarged spacing scale may be under- 
stood in terms of the (3 probability distributions. For 
a rigid rotor, with fixed /3, the spacing scale of levels 
within a band is proportional to l//3^ [1]. As observed 
above, the probability distribution of the 0^ state tends 
towards smaller than that of the ground state, giving 
a larger average 1//?^. The ratio {P~^)o+ / {P~^)o+ 
a maximum value of ^2.1 for aw 200 and decreases to 
~ 1.8 for a = 1000. Similarly, the ratio of spacing scales, 
[E{2+) - E{0+)]/E{2+), has a maximum value of - 2.3 
for a Ri 200 and decreases to 1.9 for a = 1000. 

From this understanding of the underlying mechanism, 
it is seen that enlarged energy spacing scale of the (3 band 
is an artifact of the rigid well wall in the X(5) Hamilto- 
nian [5]. As already noted, the extra energy available 
to the (3 excitation, relative to the ground state, allows 
its wave function to expand "inward" against the cen- 
trifugal potential [Fig. 5(b)], but not "outward" against 
the rigid wall. This produces the larger and hence 

rotational energy scale for the (3 excitation. A similar 
expanded spacing scale for the /? band is encountered 
in descriptions of transitional nuclei with the interact- 
ing boson model (IBM) [34] and the geometric collective 
model (GCM) [35] . In these cases the potential wall is no 



longer rigid but rather quartic (oc /3^), but the same basic 
mechanism may apply (see, e.g., Fig. 15.8 of Ref. [15]). 

Another distinctive feature of the exact solution is 
the staggering of energies within the 7 band, clearly 
visible for a ~ 200 [Fig. 2(a)]. The level energies are 
clustered as 2+(3+4+)(5"'"6+) . . ., contrary in sense to 
the (2+3+) (4+5+) . . . staggering of the rigid triaxial ro- 
tor [36]. The staggering is a remnant of the SO (5) multi- 
plet structure [28, 33] present in the 7-soft limit (a = 0) 
and is an observable manifestation of the considerable 
dynamical 7 softness seen in Fig. 4(d). It disappears 
with increasing 7 stiffness [Fig. 2(c)]. In general in a ge- 
ometric description, a weakly 7-confining potential yields 
both dynamical 7-softness and a low 7 excitation en- 
ergy, while a strongly 7-confining potential yields 7 lo- 
calization and a high 7 excitation energy. Hence, the 
presence of signatures of 7 softness (such as the 7 band 
staggering) is closely correlated with low 7 band energy. 
The quantitative relationship between these signatures 
depends upon the particular potential V{0,-y) used, so 
reproduction of the 7 band staggering may prove to be a 
valuable phenomenological test. Staggering of approx- 
imately the calculated magnitude and sense is indeed 
found in the 7 bands of rare earth transitional nuclei 
{e.g., Ref. [34, Fig. 25] or Ref. [37]). Note that resid- 
ual 7-soft staggering occurs for transitional structure in 
the interacting boson model (IBM) [38] as well, also as 
a remnant of SO (5) multiplet structure, disappearing in 
the SU(3) limit, which for infinite boson number corre- 
sponds to rigid rotor structure. 

The exact solution of the X(5) Hamiltonian has been 
seen here to provide valuable insight into the effects 
dominating /3-soft transitional structure. However, this 
Hamiltonian has several limitations as a realistic Hamil- 
tonian for detailed phenomenological analysis. 

(1) The rotational kinetic energy term in the Bohr 
Hamiltonian, cx ^^i^L'^ /[4:f3'^ sm'^{-f — ^ttk)], is con- 
structed using the irrotational flow moments of inertia, 

— ^B0^ sin^(7 — Ittk). There is extensive empirical 
evidence [1] that the actual moments of inertia are inter- 
mediate between the irrotational and rigid body values. 
This has only been established for the overall normal- 
ization of the moments of inertia, but it also calls into 
question the proper /3 dependence of these moments of 
inertia and therefore the proper (3 dependence of the ac- 
companying 7 vibrational kinetic energy term. This (3 
dependence is of central importance, since it generates 
the /3-7 coupling just found to play such a major role in 
the solution properties. 

(2) The hard wall of the X(5) square well potential 
introduces unrealistic features to the solution, as already 
discussed. The compression of the wave function against 
the well wall by the centrifugal effect was noted above 
to induce something approximating rigid (3 deformation. 
For a softer well wall, the wave function is instead free to 
expand to larger (3, so the compression should be gentler. 
The tendency of observables towards rigid rotor values for 
large 7 stiffness and the enlarged energy spacing within 
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the /? band may both be attenuated. 

(3) The restriction of V{f3, 7) to the form Vg(/3) + Ky(7) 
was imposed purely for convenience in the approximate 
separation of variables [2]. Coupling terms such as 
(3^ cos 37 arise naturally in the classical limit of the 
IBM [39, 40] and have been used extensively in earlier 
work with the geometric model {e.g., Refs. [24, 25]). The 
importance of the f3-j coupling in V{P,j) must be ex- 
plored. 

The numerical techniques of Refs. [21-23] are flexible 
and can easily be applied to a broad range of Hamiltoni- 
ans. It should therefore be straightforward to investigate 
many of the effects just discussed. Exact numerical diag- 
onalization of such arbitrary Hamiltonians may also be 
useful for more abstract studies of the geometric Hamil- 
tonian. In particular, it is well known that the geometric 
model and the IBM produce similar spectra under cer- 
tain circumstances. But much remains to be understood 
about the extent to which IBM coherent state energy sur- 
face [39, 40] can be identified with a geometrical model 
potential energy surface, and the appropriate kinetic en- 
ergy operator to be used with this potential energy sur- 



face {e.g., Refs. [41-44]). 



IV. APPROXIMATE SEPARATION 

Since extensive prior work has been carried out us- 
ing the approximate separation of variables of Ref. [2], 
it is useful to review this approximation, compare the 
results obtained under this separation with the exact 
results, and establish more clearly the reasons for the 
breakdown of the approximation. Two simplifications 
are required to obtain an approximate separation of vari- 
ables for the Hamiltonian (1). First, in the limit of 
small 7, the sum -^k / sin^(7 — 2nK/3) reduces to 
AL{L+l)/3 + K'^{l/sin^j-4/3), where L is the total 
angular momentum quantum number and K the intrinsic 
frame angular momentum projection quantum number. 
Second, in all terms involving 7, the variable (3 is re- 
placed with a constant "average" value /3o. This yields a 
Hamiltonian H = + separable in the variables /? 
and 7, 



H 



I d ,d L{L + l) 



PI 



1 d 



■ sin 3' 



d 



4 



+ V;(7). 



(3) 



The eigenvalues and eigenfunctions of Hjj are given by 
£/3 = {s = 0,1,...) and /(/3) (x r'/^ J.^/'/?) 

(0 < /3 < 1), where v = [L{L + l)/3 + 9/4]i/2 and x,,^, 
is the sth zero of Jj, [2]. Under the small 7 approxima- 
tion, the eigenproblem for reduces to that of the two- 
dimensional isotropic oscillator, with quantum numbers 
n-y and K , and the eigenvalues and eigenfunctions 77(7) 
are as described in Ref. [2]. The full eigenfunctions of 
H are products fLs{P)rin^K{'y)4>MK{'^) of a /3 ("radial") 
wave function, a 7 wave function, and a rotational or Eu- 
ler angle {■&) wave fmiction (jj^iKi"^)-: defined in terms of 
D functions as in Sec. II. 

Under these approximations, the X(5) model yields 
levels arranged in bands of good K quantum number, 
but with level energies [Fig. 9(a)] and E2 transition 
strengths [Fig. 9(b)] which differ from those of the rigid 
rotor. All predictions for the ground state band and 
/3 excitations are independent of Ky. The spin depen- 
dence of level energies within the yrast band is inter- 
mediate between those of the rigid rotor and harmonic 
oscillator, with Ri/2 = E{4^)/E{2+) « 2.90. The band 
arising from the first /3 excitation occurs at low energy 
[E{0+)/E{2t)^5.65]. 

Some of the characteristic spectroscopic features found 
in the approximate solution (discussed in detail in, e.g., 



Refs. [2, 3, 13]) are indeed encountered in the full solution 
(Sec. III). The [3 band exhibits a substantially larger en- 
ergy spacing scale than the yrast band [E{2^) — E{Op ) w 

{1.80)E{2^)]. Strong interband E2 transitions are pre- 
dicted, as is the distinctive branching pattern in which 
the spin-descending /3 to ground band transitions are sup- 
pressed. 

However, some of the essential features of the full so- 
lution are missing from the approximate solution: the 
dependence of the band head energy upon 7 stiffness, 
the confinement of the wave functions near (3=1 due 
to the five-dimensional centrifugal effect, and the con- 
sequent tendency of energy and transition strength ob- 
servables towards rotational values for large 7 stiffness. 
Also, the observable signatures of 7-softness obtained in 
the full solution for realistic 7 stiffnesses, such as the 
staggering of energies in the 7 band, are absent under 
the approximate separation of variables, which effectively 
enforces 7 rigidity in the separation process. 

The basic reason for the failure of the approximate 
solution to reproduce these features is clear from the 
analysis of Sec. HI. The approximate Hamiltonian (3) 
retains only a portion of the centrifugal (l//3^) contri- 
bution to the full Hamiltonian (1). Through the term 
L{L + l)/(3/3^), the portion of the centrifugal effect aris- 
ing from rotational kinetic energy is essentially retained. 
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Approximate 
separation 




FIG. 9: Level scheme for the approximate solution to the X(5) Hamiltonian. Excitation energies of the lowest members of 
the ground, /3, and 7 bands, normalized to E{2^), are shown at left. Electric quadrupole transition strengths, normalized 
to B{E2;2f 0^) = 100, arc shown at right. The 7 band head energy in part (a) and the 7 band transition strength 
normalizations x and y in part (b) depend upon the details of the 7 potential and are thus left unspecified. For plotting purposes, 
the 7 band head position and arrow thickness scales have been arbitrarily chosen to facilitate comparison with Fig. 2 (top). All 
observables involving the 7 band are calculated according to the approximate separation as it appears in Ref. [2] rather than 
in Ref. [3] , which would yield values differing in detail [45] . 




Radial basis flinction 

FIG. 10: Probability decompositions of the true X(5) wave functions with respect to the radial functions from the approximate 

solution, integrated over angular coordinates. Decompositions are shown for the ground, ,5, 7, and ,38 band heads, for a — 200 
(left) and a = 1000 (middle). The relevant radial wave function under the approximate separation is indicated (right) for 
comparison. A basis of L = radial wave functions from the approximate solution (i^ = 3/2) has been used for decomposition 
of the 0"*" states, while a basis of L = 2 radial wave functions {v = \fl /2) has been used for decomposition of the 2+ state. 



But the replacement of 1//?^ by l//3o in the remaining 
terms suppresses the portion of the centrifugal effect aris- 
ing from the 7 kinetic energy. 

For a more detailed understanding of the breakdown 
of the approximate separation, we can reexamine the two 



approximations made in obtaining Eqn. (3). 

(1) The small angle approximation for 7 is in principle 
arbitrarily good for sufficiently large 7 stiffnesses. How- 
ever, as discussed in Sec. Ill [Fig. 4 (right)], confinement 
to genuinely small angles requires relatively large 7 stiff- 
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nesses (a > 1000). 

(2) The validity of the other approximation, replace- 
ment of /3 with /3o in selected terms, depends upon one 
of two conditions being met. The approximation would 
be good if the wave function were sharply localized in /3, 
so indeed (}~ Pq. This condition may hold in some other 
contexts, such as small oscillations for a rigid rotor [1], 
but it is not well satisfied for the X(5) problem, which 
was specifically constructed to give a lack of (3 localiza- 
tion. Even without reference to the exact solution, it is 
seen from the approximate radial wave functions /sl(/3) 
that the the probability distribution with respect to /? is 
broad and differs significantly from state to state. For the 
approximate ground state (0^|/3~^|0f ) « 3.90, while for 
the /? band head (0^|/3"^|0_J) « 7.21, nearly a factor of 
two larger. Thus, replacing l//3^ by a "rigid" value 1//3q 
is not necessarily a small approximation. Alternatively, 
the replacement of (3 with /?o could still yield accurate 
results if the overall strength of the approximated term 
were small. This is seen from Sec. Ill to occur when the 
kinetic energy in the 7 degreee of freedom is small, and 
hence for nearly 7-soft potentials. 

The two approximations are thus valid at large 7 stiff- 
ness and small 7 stiffness, respectively. If these two 
regimes are to have any overlap, it must be at some 
intermediate 7 stiffness. Indeed, as noted above, there 
is considerable qualitative agreement between the exact 
spectrum for a = 200 [Fig. 2 (top)] and the approximate 
spectrum [Fig. 9]. However, the significance of this re- 
semblance should not be overstated, since more detailed 
calculations reveal that it arises in part from a cancel- 
lation of errors introduced by the two approximations. 
The small angle approximation tends to cause an over- 
calculation of the /? band head energy, while the rigid /3 
approximation tends to cause an undercalculaUon. 

A probability decomposition of the true wave func- 
tions with respect to the approximate radial wave func- 
tions /is, summed over angular wave functions, provides 
a direct measure of the resemblance of the true and 
approximate eigenstates. This decomposition is shown 
in Fig. 10 (left, center) for several states, for a = 200 
and 1000. The relevant radial basis state from the 
approximate solution is also indicated, for comparison 
[Fig. 10 (right)]. For intermediate 7 stiffness, the approx- 
imate solution indeed dominates the decomposition of the 
true solution, e.g., with a 96% probability for the s = \ 
basis state in the true a = 200 ground state [Fig. 10 (left)]. 
The similarity in wave functions breaks down for larger 
7 stiffness [Fig. 10 (center)]. 

In summary, the approximate separation of variables of 
Ref. [2] causes the contribution of the 7 kinetic energy to 
the five-dimensional centrifugal effect to be suppressed. 
This is always a significant contribution and is the dom- 
inant portion for large 7 stiffness. The approximate so- 
lution for the X(5) Hamiltonian qualitatively reproduces 
some aspects of /3-soft transitional structure, but only 
those which are least strongly affected by the f3-j inter- 
action. The approximate results quantitatively resemble 



the exact results only for a small range of 7 stiffnesses, in 
the vicinity of a = 200. This is approximately the 7 stiff- 
ness of phcnomcnological interest for description of the 
rare earth (iV « 90) transitional region. Consequently, 
many of the basic conclusions found in prior comparisons 
with experimental data [10-20] remain largely unaffected. 



V. CONCLUSION 

The numerical techniques of Rowe et al. [21-23] pro- 
vide a practicable approach to the exact diagonalization 
of the X(5) Hamiltonian. The wave functions and spec- 
troscopic properties obtained thereby provide insight into 
the structural features expected for /3-soft, axially stabi- 
lized transitional nuclei. The properties of the solution 
are found to be dominated by the (3-"f coupling induced 
by the kinetic energy operator, which results in a sig- 
nificant five-dimensional centrifugal effect. Most spec- 
troscopic properties are strongly dependent upon the 7 
stiffness. The results also highlight the presence of sub- 
stantial dynamical 7 softness. These basic qualitative 
features were not apparent from the approximate solu- 
tion. The analysis in the present work can readily be 
extended to Hamiltonians which provide a more realistic 
treatment of transitional nuclei. 
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APPENDIX A: MATRIX ELEMENTS 

In this appendix a summary is provided of the pro- 
cedure for calculating the necessary matrix elements of 
operators between radial or angular basis functions. The 
matrices of operators in the full product basis are ob- 
tained as the outer products of the radial and angular 
matrices. The radial basis functions /i(/3) are defined in- 
side the well (0 < /3 < 1) as fi{(3) = A^^iP-^/'^J^{x^^i0) 
and are vanishing outside the well, where Xv,i is the ith 
zero of and A„^i = [-J„-i{x„^i)J,y+i{x„^i)]~''-^'^. For 
fixed v and for i = 1, 2, . . ., these form an ortho normal 
set of functions with respect to the metric f3'^d(3. The 
angular basis functions are constructed as sums of the 
form 'i'vaLMil,^) = Y.K=o FyaLK{f)<pMKi'^) according 

even 

to the procedure of Ref. [21]. These are orthonormal 
with respect to the metric \sm3j\d'yd'&. The (ji^K ^-re 
symmetrized combinations of D functions, <pj^f({^) = 

[{2L + l)/{16n\l + 5K))mDj^Km + {-)^Di,_M*, 
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and are orthonormal under integration over the Euler an- 
gles, except that (f)j^,jj^ = in the special case oi K = 
and L odd. 

Matrix elements of an arbitrary function g{(3) between 
radial basis functions are calculated by straightforward 
numerical integration, as 

(Mi?(/3)l/i> = / P'df3fi,{(})gim{P)- (Al) 

The matrix elements of the radial kinetic energy oper- 
ator in Eqn. (1) can be reexpressed in terms of matrix 
elements of by use of the Bessel equation, as 

= xlM-{'^^-Q/^){fi'\P-'\fi)- (A2) 

Observe that for the X(5) Hamiltonian the matrix el- 
ements of V(} vanish, since Va(/3) ~ inside the well. 
The presence of the square well radial potential thus en- 
ters the calculations only implicitly, through the bound- 
ary condition it places on the allowed basis functions. 



which in turn dictates their kinetic energy matrix ele- 
ments (A2). 

The matrix elements of an arbitrary function /{"f), 
such as the function 7^ appearing in the angular poten- 
tial, can be calculated through a series of straightforward 
integrations involving the coefficients FyaLKil), as [21] 

(*^'a'L||/(7)ll*^aL) = {^v'a'LL\f{l)\^vaLL) 

= E /|sin37|rf7iC'a'LK(7)/(7)i^.aLK(7)- (A3) 

even 

Here the Wigner-Eckart normalization convention of 
Rose [46] has been used in the definition of the reduced 
matrix element. The integration can be restricted to 
0<7<7r/3by periodicity of the functions in 7 [21]. 
The expression (A3) is readily extended, by application 
of the Clebsch-Gordan series, to give the matrix ele- 
ment of any spherical tensor operator /^(7, i?), provided 
it is expanded in terms of D functions as /^(7, "!?) = 
E«=o f^{l)(l>a^{^)- The matrix element is 

even ^ 



{LKXk\L'K') 



(2£ + l)(2A+l) 
2L'-hl 



1/2 



E 

even 

' {-)^{LKXR\L'K') K>k] 
{-)^{LKXk\L'K') K<k\ 



1 + Sk' 



-,1/2 



(1 + 5k)(1+<5«) 



volLK 



(7) 



(A4) 



For the leading order E2 transition operator [1], 
mi(^2;/.)(x/3[£>^*cos7 



(A5) 



the expansion coefficients are /o(7) oc (Stt^/S)^/^ C0S7 



r 



and /|(7) oc (Stt^/S)^/^ sin7. Transition strengths are 
B{E2-i^ /) = {2Lf + l)\{^f\\m{E2)\\^i)\^/{2U + 1). 
The angular kinetic energy operator is the negative of 
the seniority operator A [28]. Since the five dimensional 
spherical harmonics are seniority eigenstates, the matrix 
elements are simply 



L'M' 



It 



1 d 



. Sin37^ Ty^^TT ^ AI^vcLm) = -v{v + 3)6v'vSa'aSL'LSM'M- (A6) 

sm37a7 07 4 sm^(7 - Ittk)/ 



Finally, consider a wave function (^(/?, 7, 1?) decom- 
posed in terms of the product basis functions as 
7, ■&) = J2i^k 0'ikfi{l3)'^kLM{l, 1?), where is a short- 
hand for the indices [va). The probability density with 



respect to /? and 7, integrated over Euler angles, is 

-P(/3,7) = /?''|sin37| ^ ^ ai'k'aik 

K i',k',i,k 
even 

X fHf3)Fk'LK{l)MP)FkLKil). (A7) 
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